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Abstract 

It is shown that if, from the starting point of a universal rank-one 
mass matrix long favoured by phenomenologists, one adds the assump- 
tion that it rotates (changes its orientation in generation space) with 
changing scale, one can reproduce, in terms of only 6 real parameters, 
all the 16 mass ratios and mixing parameters of quarks and leptons. 
Of these 16 quantities so reproduced, 10 for which data exist for direct 
comparison (i.e. the CKM elements including the CP-violating phase, 
the angles #i2j#i3,#23 in ^-oscillation, and the masses m^ra^^m^) 
agree well with experiment, mostly to within experimental errors; 4 
others {m s ,m u ,md,rn V2 ), the experimental values for which can only 
be inferred, agree reasonably well; while 2 others (m Ul ,5cp for lep- 
tons), not yet measured experimentally, remain as predictions. In 
addition, one gets as bonuses, estimates for (i) the right-handed neu- 
trino mass m UR and (ii) the strong CP angle 6 inherent in QCD. One 
notes in particular that the output value for sin 2 29\^ from the fit 
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agrees very well with recent experiments. By inputting the current 
experimental value with its error, one obtains further from the fit 2 
new testable constraints: (i) that #23 must depart from its "maximal" 
value: sin 2 26» 23 ~ 0.935 ± 0.021, (ii) that the CP-violating (Dirac) 
phase in the PMNS would be smaller than in the CKM matrix: of 
order only | sinScp] < 0.31 if not vanishing altogether. 



1 Introduction and Summary of Results 



The mass values and mixing matrices of quarks and leptons together rep- 
resent 20 independent quantities, accounting for some three-quarters of the 
total number of empirical parameters appearing in the otherwise very suc- 
cessful standard model. They fall, besides, into a quite bewildering hierar- 
chical pattern which cries out for a theoretical explanation. But while their 
experimental colleagues have now measured some of these parameters to an 
impressive accuracy and are now giving shape even to the PMNS matrix 
so incredibly difficult to access, theoreticians, embarrassingly, have failed so 
far to produce a commonly accepted explanation even for their qualitative 
features, let alone their quantitative values to the accuracy that experiments 
have achieved. 

Admittedly, the theoretician has no easy task here. The explanation for 
these parameters cannot apparently be found in the standard model itself, 
since this logically seems to admit any chosen values for these parameters. 
On the other hand, most of the current models or theories which attempt 
to go beyond the standard model tend rather to increase the number of 
independent parameters than to place restrictions on the existing ones. There 
are models giving predictions for some of the quantities, but not all. For this 
reason, the mass and mixing pattern of quarks and leptons as a whole still 
remains largely a mystery in particle physics which is still to be understood. 

The aim of this paper is to attempt a first understanding of this great 
mystery via what has become known as the R2M2 (rank-one rotating mass 
matrix) hypothesis. This starts from the assumption of a "universal" rank- 
one mass matrix for all quarks and leptons long favoured by phenomenologists 
[U E] as a first approximation and just adds to it the assumption that this 
matrix rotates (i.e. changes its orientation) with changing scale in generation 
space. A fermion mass matrix of rank one can always be rewritten by a 
suitable relabelling of the sw(2)-singlet right-handed fields [3] in the following 
form: 

m = mT<x.cx\ (1) 

and by "universal" one means that the vector ex is the same for all fermion 
species, that is, whether up-type {U) or down-type (D) quarks, or whether 
charged leptons (L) or neutrinos (N). The mass matrix (flj has only one 
massive state for each fermion species, not a bad starting point for the mass 
hierarchy, and the unit matrix for the mixing matrix (i.e. no mixing), also 
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not a bad first approximation at least for quarks, hence the original phe- 
nomeno logical interest. The rotation hypothesis in R2M2 is added to give 
the needed deviations from this first approximation, such as nonzero masses 
to the two lower generations and nontrivial mixings as observed in experi- 
ment. Notice that the fermion mass matrix should rotate with changing scale 
occurs in the standard model [3J. Thus, the only thing new in R2M2 is the 
assumption that it is the rotation which gives rise to the mass hierarchy and 
to mixing rather than the mixing which gives rise to the rotation as in the 
standard model. Practically, this means that the mass matrix has to rotate 
faster with respect to scale than it does in the standard model. 

That the mass matrix rotates means that not only the masses of the 
particles but also their state vectors in generation space have to be defined 
each at its own mass scale. That this will then lead immediately to nontrivial 
mixing and to nonzero masses for the two lower generations is easily seen 
as follows. The state vector of the heaviest state in each fermion species is 
clearly to be taken as the vector a. at its own mass scale. Thus for the t quark, 
the state vector is t = = m t ), and for the b quark, b = <*(// = m&). The 
state vectors c of c and u of u are then of necessity orthogonal to each other 
and both orthogonal to t forming together the [/-triad. Similarly, we have 
(b, s, d) the D-triad. These two triads are not aligned because the vector a 
by assumption would have rotated from the scale \i = m t to the scale \i = m b , 
hence giving immediately nontrivial mixing. Next, the vectors c and u, being 
both by definition orthogonal to t = ot(fi = mt), are both null states of the 
mass matrix ([l]) at fj, — mt- But it is not to be inferred that c and u have 
then zero masses for their masses have to be evaluated each at its own mass 
scale, not at \i — m t . Thus, the mass m c has to be evaluated at fi = m c < m t , 
at which scale, the vector cx would have rotated to a different direction and 
acquired a component in the c direction thus giving the c quark a nonzero 
mass. A similar conclusion applies to the u quark, and indeed to any quark 
or lepton of the two lower generations. The nonzero masses that they thus 
acquire by this so-called "leakage mechanism" are small if the rotation is not 
too fast, hence giving rise to mass hierarchy, qualitatively as that seen in 
experiment. 

A more detailed analysis in e.g. [5] along these lines yields the following 
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formulae for the state vectors: 



t = oc(m t )] 
c = u x t; 



u 



a(m c ) x a(mt) 



(2) 



a(m c ) x <*(m t )| ' 



and for the masses: 



m t = mu, 



in 



m. 



u 



c 



mjj\a.(m c ) ■ c| 
mu\ot{m u ) ■ u 



(3) 



Solving these equations together gives then the masses and state vectors 
of each individual state. Analogous results hold for the D-type quark and 
leptons. The mixing matrix then follows for quarks as: 



with a similar mixing matrix Uln for leptons. 

There are two points in the above result which need to be clarified, which 
the reader might have already noted. The first is that although the mass 
matrix, by assumption in R2M2, is of rank one throughout with two zero 
eigenvalues, and therefore invariant under chiral transformations in those 
eigenstates at every scale ft, yet as a result of the rotation according to (J3|, 
all quarks and leptons have acquired a nonvanishing mass. At first sight this 
may seem surprising, or perhaps even counter-intuitive, it having been in- 
grained into us by experience working with nonrotating mass matrices, that 
the eigenvalues of a mass matrix are the mass values and that chiral invari- 
ance is synonymous with the existence of zero mass states. A rotating mass 
matrix, however, of which not only the eigenvalues but also the eigenvectors 
are scale-dependent, is to most of us a new situation, and it stands to reason 
that for dealing with this we shall have to start afresh and check how many of 
those notions gained from experience before will still survive. True, the mass 
matrix is already known to rotate in the standard model 0] and the question 
of how to deal with this situation should have arisen there. However, the 
rotation there being slow, and its effects therefore negligible, practically the 
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question can be ignored, and has largely been ignored for this reason. But 
it cannot now be ignored in R2M2. Thus approaching the question anew 
with an open mind, it will become clear that for a mass matrix that rotates 
with scale, neither its eigenvalues nor its eigenvectors can in general be taken 
directly as the masses and state vectors of particles. Instead, the masses and 
state vectors are to appear as the eigenvalues and eigenvectors of a "phys- 
ical mass matrix" to be constructed in a specific manner from the rotating 
matrix. A detailed analysis of the situation is what gives the result @— Q 
above [5]. That being so, what might before have appeared as surprising or 
counter-intuitive will now no longer be the case. 

The second point to clarify concerns the Kobayashi-Maskawa CP-violating 
phase in the mixing matrix. Since the vector a is tacitly taken above to be 
reaQ all the state vectors and the mixing matrix Vud in @ deduced from 
it will be real as well. It might thus appear that R2M2 would not give 
us a proper CKM mixing matrix with a CP-violating phase. Surprisingly, 
however, R2M2 contains in itself a solution to this apparent problem, and 
in a most intriguing manner, giving at the same time a solution to the old 
strong CP problem in QCD [6]. As is well known, invariance principles al- 
low in principle in the QCD action a CP-violating term of topological origin 
with arbitrary coefficient 9, which would give rise to, for example, an electric 
dipole moment for the neutron. The non-observance of this in experiment [7j 
means, however, that this parameter 9 has to be very small, of order 10~ 9 . 
It is thus customary to insert by fiat into the standard model an implicit 
assumption that 9 = 0, or to add extra particles called axions [8] which are 
not yet observed. The need for doing so is known as the strong CP prob- 
lem. It is also well-known that the theta-angle can be removed by a chiral 
transformation on the quark fields, but at the cost of making the quark mass 
complex. Hence, unless some quark happens to have zero mass which is ap- 
parently at variance with experiment, the strong CP problem still remains. 
What is now interesting in R2M2, however, is the fact already noted in the 

2 There is in principle no need in R2M2 to assume that a is real. Even if a is complex, 
however, the relative phases between its 3 components will have to change with scale \x 
before it can give rise to CP-violations in the CKM matrix, otherwise it is equivalent to 
a real a. We have taken a to be real here because (i) we can find no empirical indication 
for the ^-dependence of its phases, (ii) we have not succeeded to build a model where a 
changes its phases with scale in the required manner, but most importantly, (iii) we see 
here that CP-violations of sufficient size will in any case arise by virtue of the theta-angle 
term in QCD even for real a. 
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above paragraph that the mass matrix can remain chirally invariant while 
keeping all quark masses nonvanishing. This means that the theta-angle term 
can be removed by a chiral transformation without making the mass matrix 
complex, while keeping all quark masses nonzero as experiment seems to de- 
mand. There is, however, a penalty, namely that the chiral transformation 
needed above to remove the theta-angle term is automatically transmitted 
by the rotation to the quark mixing matrix and makes it complex [9j. But 
this is, of course, a penalty one would be most willing to pay, for this phase is 
exactly what one was missing in the mixing matrix Vud in Q above. There 
is yet more, for it can readily be shown that provided that one starts with a 
Vud with entries roughly similar in magnitude to those of the CKM matrix 
observed in experiment, then the CP- violation transmitted by rotation from 
a theta-angle of order unity will automatically be of the same order as that 
observed in experiment, i.e. corresponding to a Jarlskog invariant [TOj [11] J 
of order 10 -5 . 

For details on the two points above and their clarification, the reader has 
to be referred elsewhere j5j [9] . Although the analysis is logically straightfor- 
ward, it does take some care and patience to sort out, for which a freedom 
from pre-conceived notions would seem, in our own experience, to be a pre- 
requisite. Some readers may find it convenient first tentatively to accept 
the R2M2 scheme as outlined above and see what it gives in practice before 
deciding to verify the details of its justification. For this reason we first sum- 
marize the results before giving the details later in Section 2 and Appendix 
A of how they have been obtained. Indeed, given the quality of the results, 
some readers may find the scheme a useful means of parametrizing the data, 
before getting into the theoretical details. 

According then to the formulae ^ — Q, one should be able to evaluate 
the masses and mixing matrices of both quarks and leptons once given ex as 
a function of the scale /i. As /x varies, the vector ex traces out a curve on the 
unit sphere, which for our purpose is best described by the Serret-Frenet- 
Darboux formulae for curves lying on a surface, as will be spelt out in the next 
section. In this language then, the trajectory for ex is completely prescribed 
first by giving the arc-length s as a function of the scale /i, and second by 
giving the geodesic curvature K g of the traced-out curve, i.e. how much it 
is bending sideways, as a function of the arc-length s. Inspired by a model 
we have been studying, we choose to parametrize s as an exponential in In \i 
(2 parameters) and K g as some Breit-Wigner shape of s (2 parameters). To 
these 4 real parameters, we have to add 6, the theta-angle from the strong 
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Quark/Lepton 


Mass used in fit (GeV) 


Experiment 


t 


172.9 


172.9 ±0.6 ±0.9 


b 


4.19 


4 1Q-HU8 
^• iy -0.06 


T 


1.777 


1.77682 ±0.00016 




0.020 


fitted parameter 



Table 1: Mass values for the heaviest generation used in the calculation 
compared with data, where for the heaviest neutrino z/3 the entry in the 
second column denotes the (unknown) Dirac mass. 

CP term to give the CKM matrix a CP-violating phase, as explained above. 
Lastly, we add the Dirac mass of the heaviest neutrino u 3 as a parameter 
for in R2M2 the rotation mechanism depends on the Dirac masses, and these 
for neutrinos are obscured by the Majorana mass term of the right-handed 
neutrino(s) via a see-saw mechanism or something similar. Altogether then, 
we have 6 real parameters to our scheme, with 9 expected to be of order 
unity. 

Once given a set of values for these parameters, the formulae ^ — Q 
then allow one to calculate all the mixing matrix elements and masses of 
both quarks and leptons in terms of the coefficients in (fi]), one for each 
fermion species, which may be identified with the mass of the heaviest state 
in that species, and given the values listed in Table [TJ From these, one 
is then to reproduce the remaining 16 quantities, namely 4 from the CKM 
matrix, 4 from the PMNS matrix (neglecting for now any possible Majo- 
rana phases), together with the 8 masses of the two lighter generations: 
m c , m s , m U2 ,m u , m d , m e , m ui . 

By adjusting the values of the above 6 parameters, as detailed in the 



next section, one arrives without much difficulty at a set (21 ) which gives the 
following results. First, for the absolute values of the CKM matrix elements, 
one obtains the following: 

0.97423 0.2255 0.00414 
|V^ M | = [ 0.2252 0.97342 0.0417 J . (5) 
0.01237 0.0400 0.999122 
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and for the Jarlskog invariant, the value: 



\J q \ = 2.91 x 10" 



(6) 



These are to be compared with experiment, for the absolute values of the 
CKM matrix elements |12|: 



0.97428 ± 0.00015 
0.2252 ± 0.0007 
0.00862l°gg§ 



0.2253 ± 0.0007 

0734^+0.00015 
u.y io^o_ 00016 

0.0403 



H0.0011 
-0.0007 



3 47+ 00016 
u - uuo ^' -0.00012 

0410 +aoon 
u.u^iu.q 0007 

9991 f;2 +0 - 000030 

u - yjy±JZ -0.000045 



(7) 



and for the Jarlskog invariant [12J: 



\J q \ = (2.911^?) x IO- 



CS) 



One notes that apart from the two corner elements, \V u b\ and \Vtd\, all other 
quantities above are reproduced to within experimental errors. The corner 
elements, being second order effects, as will be explained in the next sec- 
tion, point ( rather difficult to reproduce correctly by rotation using a 



simple parametrization such as equations (12) and (13), but have neverthe 



less roughly the right size and the right asymmetry about the diagonals, as 
already anticipated [T3] . 

For the PMNS matrix, a first question to settle is whether the chiral 
transformation applied on quarks above to remove the theta-angle term from 
the QCD action should affect leptons as well. A priori, given that leptons 
carry no colour, they need have nothing to do with the theta-angle and have 
thus no reason to undergo a similar chiral transformation. In that case, one 
would obtain for the PMNS matrix just the analogue Uln of @ above, which 
for the same values of our 6 parameters as before is found to take the form: 



th 

PMNS I 



0.818 0.556 0.147 
0.519 0.605 0.604 
0.247 0.570 0.784 



(9) 



with J e = sin 5, 



CP 



0. 



The values of the //-oscillation angles sin 2fl 12 ,sin 2 20i3, sin 2 26 23 corre- 
sponding to this PMNS matrix are listed together in Table [2] and compared 
with experiment. They are all seen to be within the experimental bounds 
given in [T2] . 



7 



Source 


sin 2 29 12 


sin 2 20 2 3 


sin 2 29 13 


Our Output 


0.864 


0.935 


0.0842 


Experiment 


0.861^;- 


> 0.92 


< 0.15 



Table 2: Output values of neutrino oscillation mixing angles compared with 
data [12] . 

Since [12] was last updated, however, several major experiments [HI [TSJ 
Uni [13 dB] have released new exciting results on # 13 which deserve a closer 
examination. This will be done in section 3 where it will be seen that our 
output value as given in Table [2] agrees very well with the new data (as 
summarized in Tables 4 and 5). 

That there is no need for a chiral transformation on leptons to remove 
the theta-angle term, however, does not necessarily mean that one cannot be 
applied nevertheless to give the PMNS matrix a CP-violating phase. This 
possibility will be examined in Appendix B where it will be shown that in 
view of the recent data on #13, the present R2M2 fit can admit only a small 
value for sin 5 l CP arising in this way. In addition, considerations will be given 
there to a situation when such a CP- violating phase may become a theoretical 
necessity. 

Lastly, the quark and lepton masses we obtained are listed together and 
compared with experiment in Table [3| The agreement is very good for c and 
/1. It is not bad even for s, e and v 2 where some delicate extrapolation is 
involved. For the others u, d and V\, no ready information is available yet for 
a meaningful comparison. A more detailed discussion on the last two points 
will be given in the next section. 

In summary, of the 16 quantities reproduced, 10 of which can be com- 
pared directly with experiment namely: 4 for the CKM matrix, the 3 mixing 
angles 6*12, #23, #13 from neutrino oscillation, and the masses m c , m M , m e , agree 
with experiment, mostly to within present experimental errors, or else to a 
high accuracy in the case of m M . The only exceptions are the corner CKM 
elements V u b, Vtd and the electron mass m e , the agreement for which is only 
approximate and outside experimental errors, for reasons that will be given 
later in the next section, points (a) and (b). Of the remaining 6, 4 can be 
compared only with values inferred from experiment, where m s ,m U2 agree 
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Quark/Lepton 


Output mass in GeV 


Experiment 


c 


1.234 


1 2 q+ u - U5 


s 


0.171* 


n inn*+ 030 

U.IUU _ .020 




0.10564 


0.10565837 




0.0064+ 


0.0083+ 


U 


0.0009* 


0.0017- 0.0033 


d 


0.0008* 


0.0041 - 0.0058 


e 


0.0007 


0.0005109989 




< 0.0016+ 


unknown 



Table 3: Output masses compared with data. Note that those entries marked 
by a * or a f cannot be compared with data directly. For the light quarks 
u, d, s marked with * the output values are supposed to represent the mass 
of the free quark measured at its own mass scale while values given in the 
PDG review [12], because of quark confinement, are determined at a scale 
of 2 GeV. For the neutrinos, marked with f, the output masses are the 
Dirac masses which are not measured directly. The "experimental values" 
cited are inferred from the measured (physical) mass differences via a see-saw 
mechanism. For more details, see text, Section 2. 
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within reason but m u , only in the rough order of magnitude, again for 
reasons to be given in the next section, points (d) and (e). This leaves 2 
others m Ul and 5 l CP for leptons as predictions, for which no experimental in- 
formation is yet available. In addition, one obtains as bonuses estimates for 
the right-handed neutrino mass m UR and the strong CP angle 9, for which no 
experimental information is likely to be available for some time, and therefore 
for the moment of theoretical interest only. Although 9 was treated above as 
an input parameter and used to fit the value of the Jarlskog invariant J q for 
the CKM matrix, that this was possible is only by virtue of R2M2 so that 
a value for 9 should itself be counted as an output. Altogether then, one 
has obtained 18 quantities for the price of 6, which is not too bad a bargain. 
What is perhaps more relevant, however, is that one has obtained, probably 
for the first time ever, a comprehensive picture of the whole mass and mixing 
pattern of both quarks and leptons with a single all-embracing mechanism 
and the same parameters throughout. 

The trajectory for ex from which these results are obtained is shown in 
Figure [Tj It is seen to possess certain characteristics, e.g. that it accelerates 
in rotation as the scale /i decreases, while the (geodesic) curvature of the 
curve it traces out is reduced, in a manner to be detailed in the next section, 
but is otherwise quite innocuous. 

More details of the calculation will be supplied in Appendix A so that, 
with the help of Figure [TJ interested readers can readily verify for themselves 
the validity of the assertions made as well as the numerical values of the 
output quantities cited. 

2 The Mechanism and How it Works 

As the scale changes, the vector a traces out a curve on the unit sphere, for 
which it is found convenient to adopt the Serret-Frenet-Darboux formalism 
for curves on a surface [19]. At every point on the curve, one sets up then 
a Darboux triad, consisting of first the normal n to the surface, then the 
tangent vector r to the curve, and thirdly the normal v to both the above, 
all normalized and orthogonal to one another. The Serret-Frenet-Darboux 
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Figure 1: The trajectory traced out by a. on the unit sphere as the scale 
/i decreases over some 6 decades in orders of magnitude from /i ~ m t to 
~ m e . The points shown on the curve are spaced at equal intervals of 1/5 
of a decade in In fi indicating from their spacing the "local" speed of rotation 
at any ji. The rotation is seen to accelerate as \i decreases, while the geodesic 
curvature (i.e. sideways bending) of the curve traced out is reduced. Also, 
along the trajectory, are marked out the individual mass scales of all the 12 
quark and lepton states to show the position they each occupy, where we 
note that u, d and e are located at the back of the sphere. 
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formulae then say: 



n' = n(s + 5s) = n(s) — k u t(s) 5s + T g u(s) Ss, 
t' = t(s + 5s) = t(s) + K n n(s) 5s + K g u(s) 5s, 
v' = u(s + 5s) = v(s) — k 9 t(s) 5s — Tgii(s) 5s, (10) 

correct to first order in 5s, a small increment in the arc-length s, where K n is 
known as the normal curvature (bending along the surface), K g the geodesic 
curvature (bending sideways) and r g the geodesic torsion (twist). For the 
special case here of a curve on the unit sphere, the surface normal n is the 
radial vector a, K n — 1 and r g = 0, so that the formulae reduce to: 

ol = ot(s + 5s) = ot(s) — r(s) 5s, 

t' = t(s + 5s) = t(s) + ot(s) 5s + K g v(s) 5s, 

v' = v(s + 5s) = is(s) — k 9 t(s) 5s. (11) 

The curve will thus be specified just by giving the geodesic curvature K g as 
a function of the arc-length s. 

For the problem here, one is interested not only in the shape of the curve 
traced out by a but also in the speed with respect to the scale /i at which 
a moves along that curve. To specify the "trajectory" of a, then, one will 
need to give not only k 9 as a function of s but also s as a function of \i. 

These two functions we propose now to parametrize as follows, inspired 
by a dynamical model [201 [2H 122] we have been considering. For s as a 
function of ln(/z) (fi in GeV) we propose an exponential, thus: 

s = exp[— a — bln(fi)] — exp[— a — bln(m t )\ (12) 

and for K g as a function of s, we propose the tail of a square-root Breit-Wigner 
function, namely: 

9 y/(s - s ) 2 + s§' ^ ^ 

with g being the residue, and so denoting the position of a pole near the 
physical region. The reason for these proposals is that in the model cited, 
there is a fixed point for ol near fi = m t , where the speed at which o rotates 
tends to zero while the geodesic curvature K g tends to infinity [22J. What 
seems to matter is only the distance of t to the pole, not the ratio of its real 
and imaginary parts, so that we can use just the one parameter sq throughout. 
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However, there is nothing very special about these parametrizations and any 
similar choice of functions would likely lead to similar results. 



With our choice of (12) and of (13) above, the trajectory for a. is specified 



by the 4 parameters a, b, g, Sq. Once given a set of values for the 4 parameters 
a,b,g,s , one can reconstruct the trajectory by integrating the equations in 



(11 ). And once one has the trajectory, one can then calculate the mass ratios 
and the mixing parameters for both quarks and leptons. 

To calculate the masses of the members of any fermion species, say the 
[/-type quarks t, c, u, one has first to normalize by supplying, say the mass 
my = m t of the heaviest state t as listed in Table [TJ The mass of the c quark 
is then obtained as per ^ and ^ by solving the equation: 

m c = mu\c ■ a(fi = m c )\ 2 . (14) 

And once we know the masses of t and c, we can then solve (|2j) for all the 
state vectors t, c, u of the [/-triad. A similar procedure applies for the D-type 
quarks to give the Z)-triad. Combining the results, we can then calculate the 
matrix Vud in @ above. 

As noted before, however, this real mixing matrix Vud is not yet the 
CKM matrix we seek. If there is a nonvanishing theta-angle term in the 
QCD action, as there is in general no reason not to have in theory, then 
strong interaction will be CP-violating, which is a situation experiment will 
find hard to accommodate. However, as was mentioned above and worked 
out in detail in [9], this unwanted theta-angle term can be eliminated in 
the present R2M2 scheme, hence making the strong sector CP-conserving 
as desired, by a judicious chiral transformation on the quark state vectors, 
thus making them complex. The easiest way to specify the effect of this is as 
follows [5]. At fi = m t , we recall that cx coincides with the state vector t for 
the t quark. The state vectors c and u must then be orthogonal to a there 
and to each other, the three quarks t, c and u being by definition independent 
quantum states. They must therefore be related to r and v, the other two 
members of the Darboux triad at /i = m t , which are also orthogonal to a, 
just by a rotation about en, say by an angle ojjj, thus: 

(t,c, u) = (a, cos uut + smujjv, cosujjjv — sinw^r). (15) 

Similarly, for the D-type quarks, we have: 

(b, s, d) = (ck', cosudt' + sina^i/, cosujdv' — suicjdt'), (16) 
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where (<*', t', f') is now the Darboux triad taken at /i = m^. What the chiral 
transformation which eliminates the theta-angle term does is to give the v 
component of each of the above vectors a phase exp(— i6/2), these becoming 
thus: 



The CKM matrix can then be evaluated with these two new triads of state 
vectors as before by Q, but this will now be the true complex CKM matrix 
containing a CP-violating phase. For instance, if one calculates the Jarlskog 
invariant^] of this new matrix, one will obtain a nonzero answer. On the other 
hand, the formulae for the masses of the various quark states are not affected 
by this chiral transformation [5]. 

The same procedure as the above for calculating the quark masses applies 
to the charged lepton states. For the neutrinos, there is a difference in 
that the R2M2 mechanism applies to the Dirac masses, which for neutrinos 
are unknown, being obscured via a possible see-saw mechanism [23] by the 
Majorana mass term for the right-handed neutrino. We have therefore to 
supply the Dirac mass of the heaviest neutrino z/ 3 as a parameter, as 
listed in Table [I] above. Once given some value for m^, the same procedure as 
before will yield the Dirac masses of the 2 lighter neutrinos. To compare these 
output values with experiment, the Dirac masses have to be inferred from 
the physical masses of the neutrinos measured in experiment by assuming 
some specific see-saw mechanism. Here, we have taken the simplest model: 
Type I quadratic see-saw mechanism [23] in which the physical masses of the 
neutrinos are given as: 



comparison with the output values, and gives incidentally also an estimate 
for the right-handed neutrino mass m UR . 

3 Note however that in R2M2, where the orientation of the mass matrix is scale de- 
pendent, CP-violation in the CKM matrix cannot be directly related to commutators 
of (hermitian) mass matrices as proposed by Jarlskog in [10 . Instead, one has to rely 
solely on the unitarity properties of the mixing matrix and work with the quartic rephas- 
ing invariants which are scale independent. These invariants appeared in earlier work on 
CP-nonconservation (11) without use being made of the mass matrix commutators. 




(17) 
(18) 
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The same procedure as the above for calculating the CKM matrix applies 
in principle also to the PMNS matrix except for the question whether the 
chiral transformation for eliminating the theta-angle term in the QCD action 
should affect the leptons as well. Here, in the text, as was explained, we 
shall consider only the case when it does not, leaving the other case to be 
considered in Appendix B. The PMNS matrix for leptons is then real and 
can be calculated in the same way as the f/D-matrix for quarks in Q above. 
As already mentioned, we neglect for now any possible Majorana phases, the 
detection of which appears to be only possible in the remote future |24j. 

Having now completed the specification of how the 16 mass ratios and 
mixing parameters are calculated given any set of values for the 6 param- 
eters a,b, g, so,9,m® 3 , we can now start varying the input values of the 6 
parameters to see whether any choice of those can give a decent description 
of the data. Before doing so, it is worth gaining first some qualitative under- 
standing of how the value of each parameter will affect the agreement of the 
output with experiment. 

Let us start with the parameters a and b, which govern the dependence 
of s on fi. Some quantities, such as the mass of the muon m M and the CKM 
element Vtb depend rather simply each on the rotation angle between two 
states, namely between r and n for m M , and between t and b for V t b- They 
would thus be sensitive to how the arc-length s, or the angle subtended 
at the centre of the sphere, depends on the scale fi, or in other words to 
the parameters a and b. The two quantities m^, Vtb have also both been 
determined accurately by experiment, and so can be used as anchors for 
these two parameters. Indeed, based more or less just on this observation, 
approximate values for a, b have already been obtained some 10 years ago 
[25J, and the new values we shall suggest will not deviate much from those. 

Next, for the two parameters g and sq governing the dependence of k 9 
on s, it is instructive first to examine the situation when 9tb, the angle be- 
tween t = ct(fj, = m t ) and b = ot(fi = m&), is small, which indeed it is 
experimentally, its value estimated from the experimental value of V t b by the 
consideration in the above paragraph being only ~ 0.041. Substituting 6 tb 



for 8s in ( 1 1 ) ,( 15 ), (16) then easily yields the following approximate formula 



for the mixing matrix Q [13J: 

cos uj — Kg sin uj 9 t b sin uj + n g cos uj 9 tb sin ujjj 9 tb 

Vud = | — sin uj — Kg cos uj 9tb cos uj — k 9 sin uj 9 t b — cosuJuO t b ) (20) 

sin to D 6 t b cos to D tb 1 
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with oj = ujo — ujjj. We notice there that the geodesic curvature n g appears 
prominently in the Cabibbo angle V us , V c d- Again, the Cabibbo angle, being 
well measured experimentally, can be used as an anchor for the parameters 



g and s appearing in (13). The same observation applies qualitatively to 
the solar neutrino angle U e2 in the leptonic equivalent matrix Uln, although 
the angles involved in this case are no longer small. Thus, by adjusting the 
parameter g, one can fit roughly the output value of the Cabibbo angle to 
its experimental value, and by adjusting the position of the pole, namely sq 



in (13), the relative magnitude of n g between the quark and lepton sectors 
and hence the relative magnitude of V us to U e2 . 

There are two things still missing in the preceding arguments, namely the 
effect of the theta-angle on the phases of the mixing matrices, and that of 
the Dirac mass of z/3. The former affects mostly the Jarlskog invariant J q for 
quarks, while the latter, by specifying where the state vector 1/3 will lie on 
the trajectory, mostly the atmospheric angle U^, and can each be used to 
anchor the corresponding quantity. However, they will both also affect the 
quantities already considered, as will also the other 4 parameters affect the 
present 2 quantities, though all to a smaller extent. Nevertheless, armed with 
these preliminary understanding of how the values of our 6 parameters will 
affect mainly which of the output quantities, it is not difficult to jiggle them 
around to achieve a tolerably decent fit of the quantities concerned. Working 
in this way, one arrives at the following set of values for our parameters, 
which is what gives the result summarized in the preceding section: 

a = 2.27635; b = 0.512; g = 0.42214; s = -0.069; 9 = 1.28; = 20 MeV. 

(21) 

There is of course nothing immutable about these parameters, nor indeed 



even about the parametrizations (12) and (13). They are meant to represent 
just a sample fit by R2M2 to data, and any parametrization which gives a 
trajectory for a with a similar shape to that shown in Figure [T] is likely to 
do quite as well. For instance, there has been no serious attempt to optimize 



the parameters in (21) on our part, for indeed, given the great disparity 
in accuracy achieved in experiment for the 16 quantities we are after, it is 
not clear what optimization of the parameters should mean in fitting them. 
What we wish to assert is that with the above sample choice, one seems able 
to obtain an adequate description of the data. 

Notice that with the 6 parameters we have permitted ourselves, we can 
anchor but 6 of our output quantities by experiment, and these 6 quantities, 
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as seen, can all be fitted without trouble to within experimental errors or to 
a high accuracy as in the case of the \i mass. But the other 10 quantities we 
also want will then have to fend for themselves. That they all have managed 
to do so reasonably well should thus be regarded as a test of the R2M2 
hypothesis. 

Before moving on further, the following notes on some points of detail 
would be in order: 

(a) First, one notes that the corner elements in (20) are proportional to 
smcuu,D Otb, where each of the angles is itself a consequence of the rotation. 
These corner elements V u b and Vu, therefore, are second order effects of the 
rotation [13] • This comes about in R2M2 just from the geometrical fact that 
for a curve on the unit sphere, the geodesic torsion vanishes so that any 
twist that these elements represent will have to be of at least the second 
order in the change in arc length 5s. It thus follows immediately that they 
are small, as experiment shows. It also follows that they have the asymmetry 
Vtd > Vub, again as experiment wants, since by virtue of the behaviour (12) of 
s as a function of /z, one has u>d > Uu- For more details on these points, see 
|13j . However, the down side is that, being second order effects in rotation, 
these corner elements are rather difficult to reproduce accurately in the R2M2 
scheme, hence their output values in (J5]) lying outside the experimental errors 
in ([7]) above. 

(b) Secondly, one notes that of the parameters in (21), 4 (a,b,g,9) are 
anchored to experiment in the high scale region with fi > 100 MeV, using 
data there of high accuracy, namely m^, V t b, V us ~ V c d and J q for quarks. 
The other 2, namely Sq and m^, are anchored on neutrino oscillation data 
of understandably lower accuracy. This means that the rotation trajectory 
is quite well determined in the region fi > 100 MeV, less well but still ade- 
quately down to about ~ 6 MeV, below which, however, it is obtained 
just by extrapolation from the high scale region. This extrapolation does not 
affect, of course, the masses of the two heavier generations, nor the two mix- 
ing matrices since these depend only on the state vectors which, including 
those of the lightest generation, are already determined by ^ at the scale 
of the middle generation. The masses of the lightest generation, however, 
depend on the extrapolation, and given the assumed exponential shape of 
(12), this can cover quite a considerable range in arc- length s, as is seen in 
Figure [TJ In view of this, that one is able to get a roughly correct value 
(i.e. 0.7 as compared to 0.5 MeV) for the e mass is already better than can 
reasonably be expected. In any case, reproducing the mass of the electron 
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being such an audacious undertaking in itself, an error of this sort probably 
needs no apology. 

(c) Thirdly, again with reference to the e mass, and those of u and d 
in the lightest generation, one notes that there are in fact more than one 
solution each to the equations ([3]). We recall that at the mass scale of the 
middle generation, say c to be specific, the state u has zero eigenvalue, but 
as the scale decreases from m c , oc(fi) will acquire a component (u| <*(//)) in 
the direction of u. This will gradually increase in size until at some point 
(u| <*(//)) matches the value \J n/m t and we have, according to (J3J) , a solu- 
tion, say the first (higher) solution, to the u mass. However, as \i decreases 
further, (u|ck( / u)) may vanish again, in which case one may find in the neigh- 
bourhood again a solution for (u\a(fj,)) = a/ fi/mt, namely a second (lower) 
solution to ([3]) for the u mass. Indeed, this is what happens to the trajectory 
here and the mass values cited in Table [3] for u, d, and e correspond to this 
lower, likely more stable solution. The limit cited in Table [3] for the Dirac 
mass of the lightest neutrino z/ 1; however, corresponds to the higher solution 
since the lower solution for has not yet been found in the range of scales 
explored and we are reluctant to trust the extrapolation further. An interest- 
ing phenomenological point to note here is that given the great disparity in 
value between the normalization factors % in Q for the 3 species U, D, L, 
namely the masses m t ,mf,,m T , one would have expected a similar disparity 
between the output masses of the lower generations obtained by leakage from 
the heaviest. This is indeed borne out for the middle generation, as seen in 
Table [3j as well as for the higher solution to u, d and e (not shown), which 
latter would be in stark contradiction to experiment. However, as seen in 
Table |3j one gets from the lower solution similar masses of order MeV for all 
3 species u, d, e, which is apparently what experiment wants. The question 
why the lightest generation, in particular, is prone to having such lower solu- 
tions to their mass in the R2M2 scheme is an interesting one deserving more 
clarification, which will be supplied later, with other details, in Appendix C. 

(d) Fourthly, quarks being confined, experiment cannot measure the masses 
of the lighter quarks s,u,d directly at their own mass scales but only at a 
higher scale, usually taken to be around 1 or 2 GeV. Since the perturba- 
tive method currently available for QCD does not allow for running to such 
low scales, there is no way at present to compare the output masses of the 
light quarks determined at their own mass scales with their masses given 
in experiment, e.g. in [12J. For the s quark, the mass scale being of order 
~ 100 MeV, not that much below the scale 1-2 GeV at which it is measured, 
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one can perhaps claim that our output value of 170 MeV at its own mass 
scale is in reasonable agreement with the value of about 100 MeV given by 
experiment at a scale of 1-2 GeV. For the u and d quarks, however, one can 
go no further than saying that the output values are roughly of the same 
MeV order as the experimentally given values. We have thus no explanation 
whatsoever so far of the crucial fact why the u quark should be lighter than 
the d. 

(e) Lastly, based on the see-saw mechanism outlined above, one can es- 
timate also the physical masses of the 3 generations of neutrinos, as respec- 
tively 0.05, 0.005, < 0.0003 eV, normalized at m U3 ~ 0.05 eV, as estimated 
from bm\ z ~ 2.43 x 10~ 3 eV 2 . Here we have normal hierarchy, as is natural 
from the leakage mechanism and universality that R2M2 implies. Notice that 
the above output values for the physical masses of v 2 and v\ would corre- 
spond to a value of 5m\ 2 of only ~ 2.5 x 10 -5 eV 2 , i.e., a factor of 3 smaller 
than the experimental value of 7.59 ± 0.21 x 10 -5 eV 2 |12j . However, since 
5m\ 2 is proportional to the 8th power of the rotation angle calculated by the 
R2M2 scheme, this factor of 3 translates only to an error of about 15 percent 
in the rotation angle, which is not bad in view of the extrapolation involved 
as detailed in (b). These masses correspond to a right-handed neutrino mass 
of about 8000 TeV, and, as far we know, there are no stringent empirical 
bounds on this [26J. 

3 A Closer Look at Lepton Mixing and the 
PMNS Matrix 

Since the PDG website [12] giving the bounds quoted in Table [2] was last 
updated, 5 experiments [MJ [15], [161 E0 HE] have given exciting new results 
on the mixing angle #13 which now need to be taken into account. We shall 
do so in the order in which these results appeared. 

As regards the two appearance experiments T2K and MINOS, the quan- 
tity they quote is 2 sin 2 6*23 sin 2 26 ) 13 , and the result depends on the CP phase 
Sq P . It is thus the quoted value of this quantity at 5 l CP = that should 
be compared with the value of 2 sin 2 #23 sin 2 2#i3 from our output above in 
Table [2} The result of this comparison is given in Table [4] where it is seen 
that our output value agrees with the experiments to the extent that the two 
experiments agree between themselves. 
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Source 


2 sin 2 e 23 sin 2 29 13 


Our Output 




0.0628 


T2K 


68%C.L. 


12 +U - Uy 


90%C.L. 


19 +alti 


MINOS 


68%C.L. 


04 +U Ub 


90%C.L. 


04 +auy 



Table 4: Output compared with the values of 2 sin 2 #23 sin 2 29 13 at b l CP = 
read from [HI E2] ■ 



Source 


sin 2 29 13 


Our Output 


0.0842 


Double Chooz 


0.086 ± 0.041(stat) ± 0.030(syst) 


Daya Bay 


0.092 ± 0.016(stat) ± 0.005(syst) 


RENO 


0.113 ± 0.013(stat) ± 0.019(syst) 


Combined Fit (Blondel) 


0.084 ±0.014 



Table 5: Output value for sin 2 2^13 compared with Double Chooz [16], Daya 
Bay [TTJ, RENO [IB] and the combined fit of the first two experiments to- 
gether with T2K [H] and MINOS [T5] by A. Blondel [27]. 

The other three disappearance experiments, Double Chooz, Daya Bay, 
and RENO measured simply sin 2 29\ 3 and the result does not depend on the 
CP phase b l CP so their result can be directly compared with our output value, 
as is done in Table [5] The agreement is seen to be very good. 

Finally one can compare our output value with the combined fit for 
sin 2 26*13 of the results from the first four experiments given by A. Blon- 
del in his summary talk at the recent Moriond meeting [27]. As seen also in 
Table [5j our output value falls right in the middle of the allowed experimental 
range. 

That our output value for #13 should agree so well with experiment is in a 
sense a little fortuitous for the following reason. As explained in the preceding 
section, of the 6 parameters for the present R2M2 fit, 5 are anchored on 
quantities which have been measured with fair accuracy, leaving only m„ , 



20 



the Dirac mass of the heaviest neutrino poorly determined, on which the 
value of #13 mostly depends. The value of does not affect much the other 
quantities under study apart from those associated with neutrinos. And of 
the elements of the PMNS matrix, it is the entries of the last column which 
it will affect the most, as can be seen as follows. The 3 elements of the last 
column may be thought of as the direction cosines subtended by the state 
vector for v% on the charged lepton triad, the squares of which add up of 
course to unity. If the trajectory of a were to lie on a plane, then 9 13 (i.e. 
9 e z) would be zero, so that the departure of 6*13 from zero can be thought 
of as a measure of the nonplanarity of the trajectory. The further the state 
vector for u 3 deviates from the plane containing ck(/i = m T ) and ot(fi = m^), 
the more sin #13 = \U e s\ = yl — |^3| 2 — | U T s\ 2 will differ from zero. Now 
we have seen already in the preceding section that in view of the largish 
values of both the Cabbibo angle \V US \ ~ \V C( i\ and the solar neutrino angle 
sin #12 ~ \U e 2\i the trajectory must be nonplanar to some fixed extent, in 
which case, the further that the vector for u 3 is from that of r, the greater 
will be the nonplanarity and the larger the value of #13 as a result. Thus, 
by adjusting the value of the parameter Tn„ 3 , namely the Dirac mass of the 
heaviest neutrino, and hence the proximity of the state vector of z/3 to that 
of r, one can change the output value of #13. It was therefore only by a bit of 
luck that by putting in ( J2l| for the value 20 MeV, a round number, that 
one has ended up with a value for sin 2 2#i 3 so close to the experimental value. 
The very good agreement obtained in Table [4] and [5] should thus be regarded 
not so much as a successful prediction, although it was indeed obtained before 
the data in Table [5] appeared, but as a test of consistency with experiment, 



first for R2M2, and secondarily for the trajectory as parametrized in (21), 
even which, of course, is already of significance. 

As noted, the output value of 9 13 can be changed by changing the input 
value of the parameter without affecting the quality of the rest of the 
fit except that as regards the neutrinos. The changes induced on the output 
values of the various neutrino quantities, however, will be correlated so that 
the output value for #13, for example, will be constrained by the experimental 



bounds on 6*12 and #23, and vice versa. Keeping then the rest of (21) the 
same and only varying m^ 3 , one readily comes to the following interesting 
conclusions: 

• sin 2 2#i3 > 0.077 if sin 2 29^ and sin 2 2#23 are to remain within their 
experimental bounds quoted in Table [2j from [12J. 
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• sin 2 2# 13 ~ 0.162 > 0.15 if #23 is to be maximal, and sin 2 2# 12 is to 
remain within its experimental bounds given in Table [2] 

In other words, both of the two conditions: (i) #13 is vanishing, and (ii) #23 
is maximal, much favoured by theoretical models previously [28] , would have 
been excluded by the present R2M2 fit, given just the earlier bounds cited in 
Table |2] from [12]. But given now the new data summarized in Tables [4] and [5] 
confirming a nonzero value for # 13 , the R2M2 fit would then give predictions 
for the amount that #23 must depart from maximal. Indeed, if one were to 
force #23 to be maximal, then the output value of sin 2 2#i3 will have to miss 
the Blondel value given in Table [5] by more than 5.5<r. Conversely, inputting 
the Blondel value for sin 2 2#i3 with his errors, then the present fit would 
predict the following value for #23 : 

sin 2 29 23 ~ 0.935 ± 0.021, (22) 

a predicted departure from maximal mixing accessible to future tests by 
experiment. 

4 Concluding Remarks 

From the results summarized in sections 1 and 3, it would seem that the 
R2M2 scheme is capable of a comprehensive description of the mass and 
mixing pattern of both quarks and leptons, which is otherwise so bewildering. 
Indeed, when supplemented by a parametrization of the rotation trajectory, 
it seems to have succeeded even in reproducing satisfactorily the mass and 
mixing parameters. To our knowledge, this is the first time that such has been 
achieved by any scheme. Given that the generations puzzle, of which these 
patterns are a part, has long been regarded as one of the great mysteries of 
particle physics, its explanation would seem a nontrivial achievement. That 
is, of course, assuming that the results and the method by which they are 
derived manage to survive the community's close scrutiny. 

It is anticipated that much of the scrutiny will be focussed on the unusual 
property of R2M2 mentioned in the introduction of being able to generate 
nonzero masses for all quarks and leptons from a universal fermion mass 
matrix of rank one. As pointed out there, this may seem counter-intuitive in 
the context of a non-rotating mass matrix and has been the cause of some 
doubts and criticisms we have encountered in private discussions. As far 
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as we are aware, these have not appeared in the public domain, although 
the R2M2 scheme has existed for more than a decade and been developed 
in numerous publications since it was proposed. But now, in view of its 
quite alluring consequences as detailed above, it would be profitable for any 
criticism to be aired and debated in public, for only when they are will the 
point be properly settled. It is after all just a matter of logic, or so at least 
it seems to us, which is capable of being sorted out. 

However, even if the R2M2 scheme manages to pass all scrutiny on its 
logical validity and all its attractive results as exhibited in this paper are 
confirmed, it still leaves the question, of course, why the fermion mass matrix 
should be of rank one, and why it should rotate. In other words, even if it is 
successful in explaining mixing and mass hierarchy, the R2M2 as it stands is 
at best an empirical rule, playing the role for quarks and leptons, perhaps, 
as the Mendeleev periodic table played for atoms or the Rydberg formula 
for the hydrogen spectrum. There is still missing the equivalent here of the 
quantum theory in atomic physics which underlies it all. 

For this, of course, one will have to go beyond the standard model, which 
leaves therefore the field wide open. We have ourselves suggested a possible 
approach based on what we call framed gauge theory (FGT) [21], namely a 
gauge theory in which are introduced as dynamical variables not only the 
usual bosonic gauge and fermionic matter fields, but also, in analogy to 
the vierbeins in gravity, the frame vectors in internal symmetry space as 
(Lorentz) scalar fields. This idea has the attractive feature of giving a ge- 
ometric meaning not only to exactly 3 fermion generations but also to the 
standard Higgs field needed to break the electroweak symmetry. It has been 
shown [2H 122] also to lead to a fermion mass matrix which is universal, of 
rank one, and rotates with changing scales, quite as needed here, although it 
is not known yet whether it will give a trajectory of the appropriate shape. 
Whatever its future, however, this framed standard model (FSM) is but one 
possible implementation of the R2M2 idea which has led to the above results. 
It is not hard to imagine many other possibilities, and some, we believe, are 
already on the way [29J. 

Optimistically then, the R2M2 scheme may play the role of, say, the 
Mendeleev periodic table to particle physics. And, as the periodic table has 
led to the quantum theory of atoms, so may, let us hope, the R2M2 scheme, 
especially now that one has gained here a rough idea of what the rotation 
trajectory should look like, give us a hint at last to solving the great particle 
physics mystery of fermion generations. 
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Postscriptum 

After this paper was completed, an article appeared on the preprint arXiv 
by G.L. Fogli et al. [32] reporting on a global analysis of neutrino mixing 
data. The values they give on the mixing angles, corresponding at la to: 
(1) sin 2 2fl 12 = 0.851 ± 0.026, (2) sin 2 26 13 = 0.096^;^, (3) sin 2 26 23 = 
0.958^qq24, are all consistent with our output values listed in Tables 2 and 
5. They note in particular the indication of a departure of 823 from its 
maximal value of 7r/4, as was predicted by our considerations in Section 3. 
The predicted value in equation (22) obtained by inputting the Blondel value 
for # 13 is within la of their fitted value. However, had we input instead 
their value of #13 from (2) above, we would have obtained the prediction 
sin 2 26*23 = 0.953^Q;g}g almost right on top of their value. They also note 
an indication for a value of b g CP = (0. 89^44) n corresponding to sin5^ P = 
0.34 ~° f 5 which is again consistent with the bound predicted in equation (27). 



Appendix A: Some computational details to facilitate 
the verification of the results obtained above. 

The results of the present R2M2 fit, as summarized in Sections 1 and 
3, are obtained from numerical calculations. In this appendix we aim to 
supply interested readers with sufficient computational details so that they 
can readily verify the results for themselves. In the numerical results below we 
provide more significant figures than necessary to avoid cumulative rounding 
errors. 

First, we give the values of the vector a at the mass scales of the various 
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quarks and lepton states: 



cr (// 


= m*) = 


(1.00000,0 


00000,0.00000) 


4- * 

<*'(// 


= rn b ) = 


(0.99912,- 


•0.04176, -0.00345) 


a 1 [fi 


= m T ) = 


(0.99762, - 


-0.06836, -0.00880) 


a 1 (/i = 


= m vz ) = 


(0.74599, - 


-0.52127, -0.41446) 




= m c ) = 


(0.99643, - 


-0.08350,-0.01281) 




= m s ) = 


(0.97052, - 


-0.22709, -0.08076) 


a (a* 


= = 


(0.95172,- 


-0.28232,-0.12048) 


a^(/i = 


= m V2 ) = 


(0.27858, - 


-0.56685, -0.77529) 




= m u ) = 


(-0.63101, 


0.76714,0.11543) 




= m d ) = 


(-0.45968, 


0.81729,0.34748) 




= m e ) = 


(-0.37336,0.81904,0.43563) 


a f (/i = 


= m Vl ) = 


(-0.71656, 


0.21263, -0.66433) 



(23) 



where we note that for consistency, the mass values of the various states are 
taken as either the values listed in Table [T] or the output values listed in 
Table pi From (23), it can readily be checked, for example by evaluating 



their polar co-ordinates, that these vectors do indeed lie on the trajectory 
shown in Figure [T] 

Secondly, we give the state vectors of the various quark and lepton states: 



tt = 


(1.00000,0.00000,0.00000) 


bt = 


(0.99912, -0.04176, -0.00345) 


rt = 


(0.99762, -0.06836, -0.00880) 


*4 = 


(0.74599, -0.52127, -0.41446) 


c t = 


(0.00000,0.98843,0.15168) 


s t = 


(0.03989,0.92267,0.38351) 


^ = 


(0.06476,0.88601,0.45913) 


u\ = 


(0.59501,0.24219,0.76636) 


ut = 


(0.00000,-0.15168,0.98843) 


dt = 


(-0.01284, -0.38331, 0.92353) 


e+ = 


(-0.02359, -0.45861, 0.88833) 


A = 


(-0.29910, -0.81831, 0.49083) 



(24) 
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It can again be readily checked that these satisfy their defining conditions in 
^ with respect to the a's in (23). 



From (24), one can then easily evaluate: 



• by taking their inner products with the corresponding en's in (23), the 
masses of the lower generation states according to the formulae given 
in @, 

• by taking the inner products between up and down states the UD 
matrix Q for quarks and the PMNS matrix for leptons. 

To evaluate the CKM matrix for quarks, one has further to supply, thirdly, 
the tangent vector r and the normal vector v to the trajectory at respectively 
fi = m t and /i = nib, thus: 

T \n = m t ) = (0.00000,1.00000,0.00000) 

i/ t (/i = m t ) = (0.00000,0.00000,1.00000) 

T t (^ = m 6 ) = (0.04179,0.98683,0.15629) 

u\n = m h ) = (-0.00313,-0.15630,0.98770) (25) 

from which one can easily check that: 

sin u v = 0.15168; sin u D = 0.23446. (26) 

This then allows one to calculate the individual CKM matrix elements, still 
just according to (HI), although, the vectors being now complex, the matrix 
elements will also be complex. The Jarlskog invariant will then be nonzero 
and can be calculated using the formula for it given in [TUl [EE] (see also 
footnote [3]) . 

In other words, with the information supplied above in this appendix, 
the interested reader can easily make spot checks on the cited numerical 
results, or even, with some patience, verify for themselves all the numerical 
results cited, and hence to satisfy themselves that the trajectory shown in 
Figure [T] does indeed lead, via the R2M2 mechanism, to the cited mass and 
mixing patterns, in close agreement with experiment, as claimed. And this, 
of course, is the primary question of our main concern. 



As to the secondary question of whether the parametrization (12) and 



(13), with the values of the parameters given in (21), does indeed give rise 



to the trajectory shown in Figure [TJ we know no other way of checking it at 
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present than by numerically integrating the Serret-Frenet-Darboux formula 
in (11), which is also straightforward to perform. 



Appendix B: The case for a PMNS matrix with a 
nonvanishing CP-violating phase (sin5 e CP ^ 0). 

As shown in [U] , because of the presence of zero modes to the mass matrix 
at every scale, the R2M2 scheme will always allow a chiral transformation 
to be made on the state in the v direction, whether of quarks or of leptons, 
without making the mass matrix complex. The effect of this transformation 
will then get transmitted by rotation to the mixing matrix and give it a CP- 
violating phase. Indeed, in the present set-up where one starts with a real 
a in (JT]), this seems the only way to make the otherwise real mixing matrix 
complex and to give the CP-violating phase a nonzero value. What distin- 
guishes quarks from leptons so far is just that such a chiral transformation 
on quarks will affect the theta-angle term in the QCD action so that this 
chiral transformation becomes now not only possible but also necessary and 
in a sense unique, given that without it, and without having the chiral angle 
matched to 9, strong interactions will be CP-violating, in contradiction to 
experiment. For leptons, on the other hand, there seems a priori no reason 
as yet why such a chiral transformation should be made, and if it is made, no 
restriction on what value the chiral angle should take. Hence, the suggested 
simplest solution in the text when no such chiral transformation is made at 
all, so the mixing matrix remains real, with sin 5q P = 0. 

There is, however, nothing in principle to stop us making a chiral trans- 
formation on the leptons and give the PMNS matrix a CP- violating phase. 
Indeed, we shall present later in this appendix even some theoretical argu- 
ments for a situation where a chiral transformation similar to that for quarks 
becomes necessary also for leptons. Let us first examine, however, what ef- 
fect such a chiral transformation will have on our previous results and what 
sort of values it can give for 5 l CP in the leptonic case. 

As was found for the quark mixing matrix earlier, one effect of the chiral 
transformation in the R2M2 scheme is to enhance the (absolute) value of the 
top-right corner element in the mixing matrix, namely, for the leptonic case 
sin #i3 . Given that this quantity is tightly constrained both by our fit and 
by the new data now available, as explained at the end of section 3, there 
will not be very much room for this enhancement to operate. Indeed, it is 
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not hard by playing with the parameters to ascertain that for # 13 to remain 
within the experimental bounds of [27] and 9 12 and #23 to remain within the 
bounds cited in Table [2] from [12], the present setup (which neglects possible 
Majorana phases) can admit only a chiral transformation on leptons with a 
value for 9 of at most 0.64, i.e. about half its previous value required for 



quarks as given in (21). This corresponds to a CP- violating phase: 

\sm8 l CP \ < 0.31. (27) 
The imposition of a chiral transformation on leptons with a 9 of the same 



value (~ 1.28) as suggested in (21) above for quarks would lead to an output 



value for sin 2 2#i3 larger than the Blondel value of Table pi by nearly 4cr. The 



conclusion (27) is of interest to experiment, but it is not easy to unravel at 
this stage how much of it is due to the R2M2 scheme and how much to the 
particular manner by which it is here parametrized. 

In the rest of this appendix we shall consider a theoretical situation where 
a chiral transformation similar to that for quarks becomes necessary also for 
leptons. This arises from the assumption in R2M2 that the vector a is 
"universal", i.e. common not only to the up and down state of quarks and 
leptons separately but common to quarks and leptons as well. Notice that 
this last assumption was not made in the original scheme of [U [2] previously 
cited, there being then no real initiative for doing so. It was made only in 
the R2M2 scheme based partly on the observation that it seems to work as 
seen above, and partly on theoretical prejudices [211 [22]. Nevertheless, once 
it is accepted, then a new situation can arise, as follows. 

The electroweak theory is usually interpreted as one in which the original 
local gauge symmetry su(2) x u(l) is broken spontaneously by the Higgs field. 
However, it could equally be interpreted, according to 't Hooft [30] and to 
Banks and Rabinovici [31], as a theory in which the local su(2) symmetry 
confines, what is broken being only a global, say su(2), symmetry associated 
or "dual" to it. In present applications of the theory, the two interpretations 
are mathematically equivalent |30| . but some may find one intuitively more 
appealing than the other. 

In the second interpretation where su(2) confines, hereafter referred to 
as the confinement picture, (left-handed) quarks and leptons, as well as the 
W — Z and Higgs bosons, appear not as fundamental fields but as bound 
(confined) states of such under su(2) confinement. For example, a left-handed 
lepton would appear as a bound state of the fundamental su{2) doublet 
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scalar field with a fundamental su{2) doublet left-handed fermion field ipL, 
forming together an su(2) singlet state as confinement requires, thus: 

Xl = 4>Hl, (28) 

while a left-handed quark would appear as a bound state of the same (p field 
but now with a fundamental su(2) doublet left-handed fermion field, say ipia, 
carrying colour a = 1,2,3, thus: 



XLa 



4>H La - (29) 



We notice, however, that the right-handed fundamental fermion fields, say 
ipR and ijjR ai a = 1,2,3, being by definition su{2) singlets, need not be so 
confined, and can thus already function as the right-handed leptons and 
quarks respectively. 

Adopting the above confinement picture for the R2M2 scenario under 
consideration, we would obtain, for the 3 basis states a,r,i^ appearing in 



(11) above, the following left-handed leptons: 

<ftat'iJ> L , 4> ] t-iI) Li $V'il> L , (30) 
and the following left-handed quarks: 

^Of^La, (t> ] T-^Lai <fi » ' *l> La, (31) 

being both bound states, under su{2) confinement, of the fundamental scalar 
with the 3-component (in generation space) fundamental fermion fields 

The question then arises where the new factors, a, r, v which distinguish 
generations, originate, i.e. whether from the bound states' scalar or fermionic 
constituent. Given that ex which appears in the mass matrix ([!]) is supposed 



in R2M2 to be "universal", i.e. the same for (30) and for ( plj ), it seems this 
factor ex. can come only from the fundamental scalar fieldy And since the 
other factors r and v are but constructs from ck, the same conclusion would 
seem to apply to them as well. 

Supposing this is true, then we have the following interesting result. We 
have noted before that a chiral transformation is to be performed on the v 



4 We note that in the framed standard model FSM [5T1[52] mentioned earlier as providing 
a possible theoretical basis for R2M2, the factor a. in the fermion mass matrix ([!]) does 
indeed originate from the scalar constituent, as stipulated. These scalar fields play in FSM 
a geometric role, namely that of framed vectors or vierbeins in the internal symmetry space. 
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component of the quark fields so as to eliminate the theta-angle term [9]. We 
note now that, in the confinement picture of 't Hooft and others as exhibited 



in (31 ), the same chiral transform can be obtained just by a phase change in 
the v component of the scalar field. That this is so can be seen as follows. 
A chiral transformation exp (20:75) applied on a quark field ip changes the 
left-handed component by the phase —a but the right-handed component by 
the opposite phase a, thus: 

exp(ia7 5 )V> = exp(za7 5 )[f(l - 75)^ + |(1 + 75)^] 

= [exp(-za)|(l - 75)^ + exp(za)§(l + 75)^] 

= exp(ia)[exp(-2ia)i(l- 75)^ + 1(1+76)^]- (32) 

Or equivalently, it changes the phase of the whole of ip by a, but the left- 
handed component by an extra phase —2a. And since a change of the overall 
phase of the quark field ip has no physical significance, this is again equivalent 
to just changing the left-handed component by a phase —2a while leaving 
the right-handed component unchanged. But in the confinement picture, we 
see that this is exactly what changing the phase of the v component of by 
—2a would do. We recall there that the left-handed quark is a bound state of 
the scalar field with a fundamental fermion, so that a change in phase in the 
former will be transmitted to the left-handed quarks. On the other hand, as 
noted above, the right-handed quarks are not bound states at all and are not 
affected by the change in phase of the scalar field. In other words, this means 
that, in the confinement picture we have adopted, just a change in phase of 
the v component of <p by —2a = 6 will eliminate the theta-angle term from 
the QCD action, and give a CP-violating phase to the CKM matrix. We 
note also that a mere change in phase of v will not affect the orthonormality 
of the basis Darboux triad. 

The really interesting thing now is the following. If the chiral transfor- 
mation to eliminate the theta-angle term is indeed to be transferred to a 
change of phase in the v component of the scalar field as above suggested, 
then the same chiral transformation would have to be performed also on the 



lepton fields, for according to (30) above, the left-handed leptons, in the con- 
finement picture, are bound states of the same scalar fields as in the quarks, 
although now to some different fundamental fermion field. As a result, the 
same argument that was applied to quarks via the rotation mechanism to de- 
duce a CP-violating phase in the CKM matrix will apply also to the leptons 
to give a CP-violating phase in the PMNS matrix. 
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However, that the "same" chiral transformation be applied to both quarks 
and leptons does not necessarily mean that the value of 9 is the same for both, 
for the two mixing matrices for quarks and leptons refer to very different 
scales and the angle 9 may itself be scale-dependent. The conclusion (27) 
reached above would thus seem to suggest, that if this theoretical scenario 
were to hold, then the strong CP angle 9 would have to lose about half its 
value in running from the quark to the leptonic scale. 



Appendix C: Solutions for the lightest generation mass 

Empirically, the lightest generation is a little anomalous. Whereas for 

the 2 heavier generations, one has m% ^> nit, > m T and m c ^> m s > m^, for 
the lightest generation one has rrid > vn u > m e and all of roughly the same 
MeV order. Of course, so long as one has no explanation at all for the mass 
spectrum, the anomaly is just a curiosity. But as soon as one aims to explain 
why the masses should take the values they do, then this departure of the 
lightest generation from the pattern of earlier generations would have to be 
accounted for. At first sight, the R2M2 scheme would seem to have difficulty 
with this, for in the mass matrix (jlj it is the heaviest generation which sets 
the scale for the masses in each species. And if the heaviest generation 
follows the above pattern, so should apparently the lower generations by 
leakage from the heaviest. This holds for the second heaviest generation, so 
why the anomaly in the lightest generation? As indicated in Section 2, point 
(c), the occurrence of multiple solutions to the mass equation ^ give us at 
least a partial answer. 

But how do these multiple solutions to the mass equation ^ come about 
for the lightest generation, and especially why only for it, not for the middle 
one as well? To see this, we recall first again how solutions to the u mass 
come about. The state vector u for u is by definition orthogonal to the tc- 
plane. At fi = m c , a.(fi) lies on the tc-plane and is therefore orthogonal 
to u and (u|o;(/i)) vanishes. As /z decreases, the vector ck(/z) rotates away 
from the tc-plane giving then a nonvanishing component (u|ai(yu)) in the u 
direction, which gradually increases until at some value of fi it matches the 
value \J iijvn% and one has a solution to the equation (J3J) for the mass of 
u. This is then our first solution to the u mass, which would occur on the 
trajectory not so far from c, and has in actual fact a value of about 300MeV, 
which is quite different from what experiment indicates. 
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However, as fi lowers further, especially when the arc length s is growing 
exponentially with fi as suggested by (12), it can easily happen that the 
trajectory should hit the tc-plane again. Then at this point, one has again 
(u|o;( / u)) vanishing, and the same opportunity as before will present itself 
for again a solution for the u mass. This is exactly what happens for the 
trajectory of Figure [T] and accounts for the output mass value for u listed in 
Table |3j One sees therefore that it is rather easy for the u mass to acquire 
this lower solution; all it needs is for the trajectory to hit the tc-plane again 
at some lower scale /i, which it can easily do so long as the trajectory is not 
prematurely terminated, say by a rotational fixed point, for example, before 
it has gone that far. 

The same argument would hold also for d and e if one just replaces the tc- 
plane above by respectively the bs- and r^-plane. These planes are defined 
at the higher scales of the two heavier generations where the trajectory is 
moving relatively slowly so that the planes are not too different from each 
other in terms of the arc length s, nor from the tc-plane for u. Hence the 
lower solutions for u, d, e would occur consecutively and fairly close to one 
another on the trajectory, as is indeed seen in Figure [TJ And since at these 
lower scales, the trajectory is moving fast, meaning that changing s requires 
very little change in fi, the solutions for the masses of u, d, e would be rather 
close in values to one another, as is indeed seen in Table [3] So much then 
for the properties noted in the text for the output masses for the lightest 
generation. 

A natural question to ask at this point is, of course, what about the 
middle generation; will there not be other lower solutions to their masses 
too? Interestingly, the answer is no, there being a rather subtle difference 
in the R2M2 scheme between the 2 lighter generations which makes the 
likelihood of lower solution practically zero for the middle generation. This 
comes about as follows. We recall, as we have done already several times, 
that all the state vectors are determined in the R2M2 scheme by the scale 
of the middle generation. This means in particular that by the u mass scale 
where one looks for solution to the mass equation (|3| for m u , the state vector 
u is already known, namely in that direction orthogonal to the tc-plane. In 
contrast, in seeking a solution to the c mass at scale fi = m c , one has first to 
ascertain what is the state vector c lying on the plane orthogonal to the state 
vector t of t. In words, what the equations ^ and ^ say is that one should 
project the vector ot(fi) at fi = m c on to the plane orthogonal to t. Then 
the direction of the shadow cast by cx(m c ) on to that plane is the vector c, 
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and the length of that shadow squared times m t is the mass m c . In other 
words, for a solution to the c mass at a scale fi, one will need the shadow 
to be of order only y/jjjmt, a small number especially when jj itself is small 
when looking for a lower solution. This means that for a lower solution, one 
will need to have <*(//) nearly parallel or antiparallel to t itself. In other 
words, for m c to have a lower solution, one will need the trajectory to pass 
again at a lower scale very close to the starting point t = (1, 0, 0) or to its 
antipodes (—1,0,0), which can happen, of course, but is very unlikely, and 
it did not happen for our special trajectory in Figure [Tj The demand is very 
much more stringent than for a lower solution for the u mass which needs 
only that the trajectory hits again a plane, i.e. the tc-plane, which it will 
almost always do so long as it is not terminated prematurely. The difference 
comes about just because there are 3 generations in nature and generation 
space is 3-dimensional. 

We note that the trajectory in Figure [T] actually passes through the tc- 
plane at fi around 1 MeV so that there are two solutions very close to each 
other to the mass equation ^ for the u mass, one lying just above and 
one just below the tc-plane. Numerically, they differ by only about 2.5 
KeV, and to the accuracy one is interested in at present, may be regarded 
practically as the same solution. Theoretically, however, we are yet unclear 
what this double solution means. We note also that if we further extrapolate 
our trajectory, it might pass through the tc-plane again and give further 
and even lower solutions to ^ for the u mass. We are reluctant to do so, 
however, for we do not believe that the exponential growth in rotation speed 
as parametrized by (12) can go on forever. We think it is more likely that 
the trajectory will eventually terminate at some low scale fixed point as it 
actually does in the theoretical model for R2M2 [2"Tj |2"2"] that we are studying. 
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